AD-A075  455 


PRINCETON  UNIV  N  J  DEPT  OF  CIVIL  AND  6E0L06ICAL  ENG— ETC  F/G  8/12 
THE  CRITICAL  VELOCITIES  OF  A  FLOAT I N6  ICE  PLATE  SUBJECTED  TO  IN— ETC (U) 
AUG  79  A  D  KERR  DACA89-73-C-0228 


NL 


UNCLASSIFIED 


CRREL-79-19 


REPORT  79-19 


The  critical  velocities  of  a  v— 
floating  ice  plate  subjected  to 
in-plane  forces  and  a  moving  load 


CRREL  Report  79-19 


V/ 


The  critical  velocities  of  a  floating  ice  plate 
subjected  to  in-plane  forces  and  a  moving  load 


I 

Arnold  D.  Kerr 


August  1979 


i 


r 

I 

I 

r 

( 


Prepared  for 

DIRECTORATE  OF  MILITARY  PROCRAMS 
OFFICE.  CHIEF  OF  ENCINEERS 

#Y 

UNITED  STATES  ARMY 
CORPS  OF  ENCINEERS 

COLD  RECIONS  RESEARCH  AND  ENGINEERING  LABORATORY 
HANOVER.  NEW  HAMPSHIRE,  U  S  A. 


Approved  for  puMk  release,  distribution  unlimited. 


Unclassified 


SECURITY  CLASSIPICATION  OP  THIS  PAOB  (Whmt  M<  BM-rmO 


REPORT  DOCUMENTATION  PAGE 


SlaMIM-) 


THE  CRITICAL  VELOCITIES  OF  A  FLOATING  ICE  PLATE 
^bJIcTED TOmPLANE  FORCESAND  A  M<5VIn!£lOAD  4 


.  AUTHOR!*} 

Arnold  Xk" 


NIZATIOM  NAME  AND  ADDRESS 

U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory 
Hanover,  New  Hampshire  037SS 


n.  control  li  no  or  net  name  and  address 
Directorate  of  Military  Programs 
Office.  Chief  of  Engineers 
Washington,  DC  20314 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


S.  RECIPIENT'S  CATALOO  HUMSER 


S.  TYPE  OP  REPORT  A  PERIOD  COVERED 


a.  PERPORMINO  ORO.  REPORT  NUMBER 


s.  contract  or  orant  NUMBER!* 


Grant  DACA  89-734)228 


NCY  NAME  4  ADDRESS!!!  dllloiant  from  Controlling  Olllc a)  I  IS.  SECURITY  CLASS,  fol  Mils  topot i) 


Unclassified 


a.  distribution  statement  (oi « m*  sipvo 


Approved  for  public  release; distribution  unlimited. 


W 


QACA8 9-7 3-C -0228 


•7.  DISTRIBUTION  STATEMENT  (ol  M<  aA*ir-cl  «Mna  In  Black  so,  II  dlltaronl  tram  Baporl) 


IS.  KEY  POROS  (CanllMM  an  n—f«<  iM  II  nKiaa—T  ana  Igantlty  By  block  numbar) 

Critical  velocity 
Floating  ice 
Ice 

Loads  (forces) 


ABSTRACT  fCanMma  an  raaocoo  alR  l(  nacaaaary  mg  Iganllfy  by  block  mmbot) 

The  critical  velocities  of  loads  moving  over  floating  ice  plates  have  been  determined  by  several  authors.  In  all  these 
analytes  it  was  assumed  that  the  in-plane  force  field  in  the  Ice  cover  is  zero.  However,  due  to  constrained  thermal 
strains,  in-plane  forces  do  occur  in  the  field.  The  purpose  of  the  present  paper  is  to  determine  their  effect  upon  the 
critical  velocities  of  the  moving  loads.  It  is  shown  that  a  uniform  compression  force  field  reduces  the  critical 
velocity,  whereas  a  tension  force  field  has  the  opposite  effect. 


I  jSTts  1473  edition  or  I  NOV  ••  is  oosolete  ^  £  Unclassified 


SECURITY  CL ASSIPIC AVION  OP  THIS  PAGE 


ii 


^ f  [  i  I  ’  PREFACE 

This  report  was  prepared  by  Dr.  Arnold  D.  Kerr,  Visiting  Professor,  Department  of  Civil  Engineer¬ 
ing,  Princeton  University,  under  Grant  DACA  89-73-0228  from  the  U.S.  Army  Cold  Regions  Re¬ 
search  and  Engineering  Laboratory.  The  work  was  performed  under  DA  Project  4A062 103 A894, 
Engineering  in  Cold  Environments.  Task  02,  Engineering  Design  Criteria,  Work  Unit  005,  Crossing 
of  Water  Barriers. 

Technical  review  of  this  report  was  performed  by  Guenther  Frankenstein  and  Dr.  Donald  Nevel  of 
CRREL. 


THE  CRITICAL  VELOCITIES  OP  A  PLOATING  ICE  PLATE 
SUBJECTED  TO  IN-PLANE  FORCES  AND  A  MOVING  LOAD 

by 

Arnold  D.  Kerr 


INTRODUCTION 

It  Is  well  known  that  when  a  vehicle  is  movint  on  a  floating  ice  sheet  there  is  a  certain 
velocity,  denoted  in  the  mechanics  literature  as  the  "critical  velocity"  vcr,  at  which  it  may  break 
through  the  ice.  Methods  for  determining  vct  for  floating  plates  subjected  to  moving  loads  have 
been  presented  by  Assur  (1961).  Kheishin  (1963.  1967)  and  Nevel  (1970).  Kheishtn  and  Nevel 
utilised  the  linear  bending  theory  of  plates  to  describe  the  response  of  the  ice  cover,  and  the  equa¬ 
tions  of  an  ideal  fluid  to  describe  the  response  of  the  liquid  base. 

Recently.  Kerr  (1972)  showed  that  an  axial  force  in  a  beam  on  a  Winkler  base  that  is  subjected 
to  a  moving  lateral  load  may  have  a  profound  effect  upon  vcr.  Since  axial  in-plane  forces  also  occur 
in  floating  ice  covers,  caused  by  constrained  thermal  strains,  it  is  of  interest  to  determine  their 
effect  upon  vcr.  The  purpose  of  the  present  paper  is  to  study  this  phenomenon. 

At  first,  we  analyze  some  related  problems:  the  magnitude  ot  the  critical  in-plane  compression  forces 
for  a  uniform  biaxial  stress  field,  and  then  the  propagation  of  free  waves  in  the  floating  plate  subjected  to  in* 
plane  stresses.  This  is  followed  by  the  determination  of  vC(  for  a  floating  plate  subjected  to  a  moving  load 
and  an  in-plane  force  field. 


THE  INSTABILITY  OP  A  FLOATING  INFINITE  PLATE 
SUBJECTED  TO  A  BIAXIAL  FORCE  FIELD 


It  is  assumed  that  the  governing  equation  for  the  determination  of  Nct  is 
DV*w  ♦  NVew  +  pjgw  =  0 

where 


V4 


At4  AxW  <9y4 


(1) 


(8) 


w(x.  y)  is  the  perturbation  from  the  plane  state.  D  is  the  flexural  rigidity  of  the  plate,  and  Pfg  is 
the  specific  weight  of  the  liquid  base. 
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Figure  1. 


For  a  buckling  mode  of  the  form 


w(x.y)  -  w„  sin(a,x)  (3) 

where  wQ  const.,  eq  1  yields,  setting 
P($/D  =  k4  =  l//‘. 

("»  '  ^  al  +  K<)wo  =  0.  (4) 

Equation  4  is  satisfied  for  a  non-zero  w  when 

«?  “  £  al  ♦  *4  =  0  (5) 


and  thus  when 


N 

D 


(6) 


Since  the  right-hand  side  of  eq  6  is  N  0  it  follows,  as  anticipated,  that  the  deformation  mode  of 
the  form  shown  in  eq  3  is  only  possible  when  N  is  a  compression  field,  as  is  shown  in  Figure  1. 
From  b(N/D)l&a}  =  0.  it  follows  that  N  =  when 


a,  =  «. 


Substituting  eq  7  into  eq  6.  it  follows  that 


cr 


N 

—  2  k*  =2  W  — 

D  \l  D 

or,  rewritten, 

N„-ly/&BT 

The  corresponding  wave  length  is 
2rr  2e 


\  -  —  -  —  -  2 nil  =  2* 


*i 


(7) 


181 


(9) 


(10) 


THE  PROPAGATION  OP  FREE  WAVES 


The  liquid  layer 

Assuming  that  the  liquid  layer  (Fig.  2)  responds  like  an  ideal  fluid  (Kheishin  1963),  and  that 
the  ratio  of  amplitude  to  wave  length  is  very  small,  the  resulting  equations,  in  terms  of  the  velocity 
potential  9.  are: 


*  Th«  urn*  result  I*  obtained  tor  w(x)  «  **>„  tin  (o,x)  tin  However  then,  instead  »f  e«i  T,  s  j  *  a  j  ♦  a j. 
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aurfOCt 


where 


Assuming 


and  substituting  it  into  eq  11.  it  follows  that  <f>(z)  has  to  satisfy 


The  general  solution  of  the  above  differential  equation  is 


Noting  boundary  condition  13.  it  follows  that 


4kz)  =  2 BleaH  cosh  a  (H  -  z). 

Substitution  of  eq  15  and  14  into  the  remaining  boundary  condition  (12)  yields 


The  above  equation  is  satisfied,  tor  any  amplitude  of  when 


tgh(aH) 


Noting  that  a  =  2n/\  it  follows  that  a  solution  of  the  form  14  will  exist  when  the  phase  velocity  c 
and  the  wave  length  A  are  related  by  the  dispersion  relation: 


(16a) 
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Equation  16  ia  presented  schematically  in 
Figure  3.  Thus,  there  exists  a  range  of  0  s  c  <  ~ 
and  to  each  phase  velocity  there  corresponds  a 
wave  length  A. 

It  may  he  shown  (Lamb  1945)  that  to  the  above 
solution  there  corresponds  a  sinusoidal  (tee  surface 
wave  of  wave  length  A  which  travels  in  the  x 
direction  with  the  velocity  c.  For  additional  com¬ 
ments  the  reader  is  referred  to  the  literature  on  (tee 
waves  in  liquids. 


The  stressed  plate  (withoet  base) 

Assuming  that  the  plate  responds  elastically 
and  is  subjected  to  a  uniform  in-plane  force  field 
N.  the  following  plate  equation  is  used  for  long 
waves: 


Assuming  a  traveling  wave  of  the  form 


and  substituting  it  into  the  above  differential  equation,  it  follows  that 
(Da4  -  Wa8  -  Pphaec8lw0  sin  a  (x  -  ct)  »  0. 

Thus,  s  wave  of  the  form  1 8  can  propagate  in  the  plate  when  N  <  Da2  and  then  c,  a  and  N  are  related  by  the 
dispersion  relation 


Writing  W  ■  na2D  where  n  «  I ,  eq 


Equation  19a  is  presented  schematically  in  Figure  4. 


Tbs  atwaaed  plats  renting  ea  a  liquid  layer 

Retaining  the  assumptions  made  in  the  two  preceding  sections,  the  resulting  formulation  is 
(Kheiahin  1953): 


.  ... 
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(21) 

<*> 

dw 

to 

*=W  fit 

(22) 

d* 

to 

=  0 

s-H 

(23) 

and  regularity  conditions  for  w  and  ♦  as  x*  +  ye  -  » . 

It  is  assumed,  as  above,  that 

w(x.y.  /)  =  w0  sin  a  (x  -  cl)  (18) 

♦(*.  I,  t)  =  A<f,(t)  cosa  (x  -  ct).  (14) 

Note  that  to  the  assumed  4>  in  eq  14  there  corresponds  in  the  liquid,  and  hence  also  at  the  piate/liquid  inter¬ 
face,  a  sinusoidal  wave  of  the  type  shown  in  eq  1 8. 

It  was  shown  before  that  differential  equation  81  and  the  boundary  condition  83  are  satisfied 
when 

<£(*)  =  28,  eaH  cosho  (H  -  x).  (15) 

Substituting  expressions  IS  and  18  into  condition  22  yields 


i-SABje^Isinhatff  -  x)]#=w  ♦  cw0l  a  cosa  (x  -  ct)  =  0. 


The  above  equation  is  satisfied  when 


2ABl 


cv0 

e®*  sinha  (If  -  w)  ‘ 
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Thus 


cw0  cosh  a  (H  -  x)  cw„  cosh  a  (H  -  x) 

♦(* .*.*)  -  —  - - - cosa  (*  -  cO  a - - — - cosa  (x  -  ct)  (24) 

sinha(H-v)  sinha  H 

since  usually  w  <<H. 

Substitution  of  w(x.t)  as  given  in  eq  18  and  of  ♦(x.x.t)  as  shown  above  into  differential  equa¬ 
tion  20  yields,  noting  that  w  «H, 

l Da4  -  Nas  -  i>  hate*  ♦  Pli  -  p{  a~  ^°S^a  W)s>0  sin«  (x  -  ct)  -  0. 

'  y  sinn  a  H  ' 

The  above  equation  is  satisfied  when 

a1a8-^fTk;cl,*H  0  l85) 

This  is  the  dispersion  relation  for  the  stressed  plate  on  a  liquid  base. 

Thus. 


(Da8  -  N)  ♦ 
_ 

lP  atghtaH)  | 


(26) 


As  expected,  when  the  plate  is  absent  (D  =  N  =  h  -  0)  eq  26  reduces  to  eq  16.  and  when  the  liquid 
is  absent  (pf  -  0)  eq  26  reduces  to  eq  19. 

Denoting 


eq  25  may  be  written  as 

1  ♦  (at)4  -  — -  ^pf  +  atctgh  (a^)|  •  0 


(25a) 


and  eq  26  as 


vj  ■  I ♦<«<)* M)2  -M'm  ,26a) 

*  (o«)  \+nattgh(aH)  ' 

A  graphical  plantation  of  eq  26  it  shown  in  Figure  S  for  H  *  «>,  pt  *  I  g/cm,,pp  ■  0.92  g/cm*,  F  * 

SO ,000  kg/cm1,  v  ■  0.34  and  h  *  30  cm  (thus,  8  *  4S4  cm  and  n  ■  0.0609).  The  situation  is  similar  to  that 
of  a  beam  on  a  Winkler  base  discussed  recently  by  Kerr  ( 1972).  Namely  for  a  given  N  <  Ntt  a  wave  of  the  form 
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18  cannot  propagate  fot  0  <  c  <  (cpf)mta  and  foi  each  c  >  two  wave  trains  with  different  X  may 

propagate. 

The  findings  of  Kerr  (1972)  as  well  as  those  of  Kheishin  (1967)  suggest  that  (c  f)m|  may  be 
the  critical  velocity  of  a  load  P  which  moves  in  the  direction  of  the  x-axis.  m  " 

A  necessary  condition  for  locating  (c  .)  i  ,  for  a  given  N  is 


It  yields,  using  eq  26.  and  then  setting  \/a*D  *  n  to  simplify  the  presentation  of  the  final  result 
ll  *  («Z)4(1  -  n)]aH  -  1 1  -  (o/>4|  smh3(o//)+  H  sinh  (2ciW)I(«/>4(3  -  n) 


or.  rewritten 


1  >  (ol>4(ff  -  3)  2aH  [l  -  (ab^Z^abtghiaH) 
1  ♦  («/)4(l  -  n)  si"h(2a«)  1  +  («j£)4(l  -  n) 


It  is  of  interest  to  note  that  eq  26a  and  27a  fa  the  case  N  0  are  the  conditions  fa  the  determina¬ 
tion  of  the  critical  velocity  vCf  used  by  Nevel  (1970). 

The  minimum  value  of  N  is  obtained  (torn  the  condition 
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It  is  found  to  uke  place  at 

a  « 

The  corresponding 

WmiJcpj-0  B  2  >/?,&.  (28) 

Thus,  it  is  equal  to  the  critical  ccmpiession  force  Na  given  in  eq  9. 


THE  STRESSED  FLOATING  PLATE  SUBJECTED  TO  A  MOVING  LOAD  P 

It  is  assumed  that  the  plate  is  subjected  to  a  force  P  which  moves  with  a  constant  velocity 
Vq  as  shown  in  Figure  6.  The  governing  equations  are: 


PCj  v*  ♦  NVl  yw  «•  p  h  —  +Pjgw  -  pf  =  P$(x  -  v0/)5(y) 

y  ,y  v  *2  <*  *=» 


-o 


(29) 


where 


dw 


dz  la=w  (N 


(30) 


9*\ 

—  *  0 
3z  la=H 

and  regularity  conditions  for  w  and  <&  as  x8  +  y2  -  «• . 


(31) 
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it  is  further  assumed  that  after  a  time  t  a  steady  state  will  exist.  This  assumption  allows  the 
time  variable  to  be  transformed  out  by  means  of  the  simple  transformation 

£  =  *  -  v  f  =*  y  <=z- 

The  resulting  formulation  consists  of  differential  equations 


with  the  boundary  conditions 


(94>  <9w 

£=w  V°  <9£ 


<9d> 

K 


<=H 


=  0 


and  the  regularity  conditions  for  w  and  d>  as  x2  +  y2  -  ~ 


lim  \  .0. 

£>7-»  +  <»  V  <9£  «9?y  <9£  <9q  / 


(32) 

(33) 


(34) 

(35) 


(36) 


The  resulting  formulation  differs  from  the  one  used  by  Kheishin  (1963,  1967)  and  Nevel  (1970) 
in  that  it  contains  the  term  NV*tri  w  in  the  first  equation.  To  solve  the  posed  problem  we  proceed 
according  to  Kheishin  (1967)  and  Nevel  (1970). 

We  introduce  the  double  Fourier  transforms 


w(a./3)=  J  j  d£dv  (37) 

—oe  —00 

*(a.(3,0=J  f  Wt.V.O^+^dtdr,.  (38) 


Multiplying  eq  33  by  e‘,(a^+  and  integrating  it  as  indicated  in  eq  38  we  obtain 


7  r wo- 

_»  _»  y<9£®  <9»/2  J 


Using  integration  by  parts,  and  noting  regularity  conditions  36 
tion  becomes 


and  definition  38.  the  above  equa¬ 


tor  0  S  <<  ~ 


v:‘ *  -I*  ^ 


^-^*  =  0 


(39) 
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where 


y*-  a* 

The  same  transformation  performed  on  the  boundary  conditions  34  and  35  yields 


<K 

<*> 


■  iav0W 


C» 


0. 


The  general  solution  of  eq  39  is 


«1>  + 


Noting  the  boundary  conditions  in  eq  41  and  42.  ♦  becomes 

*( C08hy(W  ~  0 
y  sinh  y{H  -  w) 


For  the  aliove  equation  simplifies  to 

<Ma.  ft.  £)  a 


-Ol 


y  sin  yH 

The  same  transformation  performed  on  differential  equation  32  yields 

|D(u4  4  2<i*ftP  4  ft4)  -  N(„*  4  ft*)  -  phvgn8  4  \w  +  if»fv0nd*  I  P. 

I  V  " 


140) 

(41) 

(42) 

(43) 

(44) 

(45) 

(46) 


Since  according  to  eq  45.  for  w  ^  H, 
iav,, 


d» 


ytghyH 
eq  46  becomes,  noting  eq  40. 


0  v 


Dy4  -  Ny*  4  ,t(g 


vn  (ptfln*  ♦  — — — )  W 
0  V  ><thyH/\ 


P, 


Denoting 


t4^ 

Pt* 


Pvh 


v8 

L41!  V0 

0  it 


(47) 


eq  47  may  be  rewritten  as  follows 


— 
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P 


w 


Pf 


Inverting  w.  it  follows  that 


Mf.  n) 


P  7  f  +  Pi* 


ow  40 

If 


where 


4n^pfg  ba^/y 


i  1  +  ( yl )*  -  — —  (yhft 


da  dp 


vl  +  yt*&(rtj)\  • 


u 


It  may  be  shown,  as  done  by  Nevel  (1970),  that  w({,  rj)  -*  «*>  when 
a  *  b 


and 


(48) 


(49) 

(50) 

(51) 

(52) 


Nevel  used  these  two  conditions  to  determine  the  critical  velocity  v0. 

It  Jiould  be  noted  that  these  two  conditions  are  identical  to  eq  25a  and  eq  27a,  when  the  variable  (yC)  is 
replaced  by  (<rt)  and  the  parameter  V0  by  Vpf.  Since  the  variables  (yS)  or  (a8),  determined  from  eq  52  or  27a, 
are  identical,  it  follows  that  eq  52  for  the  determination  of  V0  and  eq  25a  for  the  determination  of  Vpf  are 
identical.  Thus,  as  expected,  the  critical  velocity  v„  is  the  same  as  (cpf)mh>.  Hence,  the  dependence  of  the 
critical  velocities  v„  upon  the  axial  force  field  N  is  as  shown  in  Figure  5  for  (cpf)min  and  is  shown  in  Figure  8 
for  h  *  30  cm  and  90  cm.  The  dependence  of  vCT  upon  N  and  h  was  obtained  by  numerically  evaluating  eq  26 
and  27  for  //**••,  noting  that  (cpf)mto  =  va.  The  obtained  results  are  shown  in  Figure  9. 

— 
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CONCLUSION 

The  effect  of  a  uniform  in-plane  force  field  in  a  floating  plate  upon  the  critical  velocities  of  a 
moving  load  has  been  studied.  For  an  increasing  compression  force  field  the  critical  velocity  de¬ 
creases,  approaching  the  value  zero  as  N  -  NCf  ;  tot  an  increasing  tension  force  field  vcr  increases. 
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